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Abstract 

Fewnomial theory began with exphcit bounds — solely in terms of the number of 
variables and monomial terms — on the number of real roots of systems of polyno- 
mial equations. Here we take the next logical step of investigating the corresponding 
existence problem: Let FEASjk denote the problem of deciding whether a given sys- 
tem of multivariate polynomial equations with integer coefficients has a real root or 
not. We describe a phase-transition for when m is large enough to make FEASk be 
NP-hard, when restricted to inputs consisting of a single n-variate polynomial with 
exactly m monomial terms: polynomial-time for m<n + 2 (for any fixed n) and NP- 
hardness for m>n + n^ (for n varying and any fixed e > 0). Because of important 
connections between FEASk and ^-discriminants, we then study some new families of 
j4-discriminants whose signs can be decided within polynomial-time. (j4-discriminants 
contain all known resultants as special cases, and the latter objects are central in algo- 
rithmic algebraic geometry.) Baker's Theorem from diophantine approximation arises 
as a key tool. Along the way, we also derive new quantitative bounds on the real zero 
sets of n-variate {n + 2)-nomials. 



1 Introduction and Main Results 

Let FEASk — a.k.a. the real feasibility problem — denote the problem of deciding 
whether a given system of polynomial equations with integer coefficients has a real root or 
not. While FEASr is arguably the most fundamental problem of real algebraic geometry, 
our current knowledge of its computational complexity is surprisingly coarse, especially for 
sparse polynomials. This is a pity, for in addition to numerous practical applications (see, 
e.g., [BGV03J), FEASk is also an important motivation behind effectivity estimates for the 
Real NuUstellensatz (e.g., |Ste74| ISchOOj ) . the quantitative study of sums of squares |Blen4j . 
and their connection to semi-definite programming |Parn8j . Furthermore, efficient algorithms 
for FEASk are crucial for the tractability of harder problems such as quantifier elimination, 
and computing the closure and frontier of more general types of varieties such as sub-Pfaffian 
sets iGV04. 
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Before stating our main results, we will need to clarify some geometric notions concerning 
sparse polynomials. 

Definition 1 Let f{x) := ^f=iCjX"-^ G Z[xf^, . . . ,x^^], where x"-^ := x^^^ ■ ■ ■ Xn"\ Cj ^ 
for all j, and the aj are distinct. We call such an f an n-variate m-nomial and we 
call {ai, . . . , ajn} the support of f . Also, for any collection T of polynomial systems with 
integer coefficients, let FEASk(^) denote the natural restriction o/FEASr to inputs in T . o 

Note that n-variate quadratic polynomials are a special case of n-variate m-nomials with 
m = 0(n^). 

Definition 2 Let us use for the cardinality of a set S and say that a subset A = {ai, . . . , 0^} cZ" 
with ^A = m is afiinely independent iff the implication '\YlY=i1j'^j — O SjLi7i~0] 
=^ 7i= ■ ■ ■ = 7m = 0" holds for all (71, . . . ,7m) gM™. Also let dim A denote the dimension 
of the subspace o/M" generated by the set of all differences of vectors in A. o 

Note in particular that i^A> 1 -|- dim A, with equality iff A is affinely independent. It is also 
easily checked that a random set of m points in R" (chosen, say, independently from any 
continuous probability distribution) will have dimension min{n, m — 1} with probability 1. 

Clearly, any n-variate Laurent polynomial / always satisfies dim Supp(/) < n,. Let M*: = 
M \ {0}. It is not much harder to see that one can always find (even in an algorithmically 
efficient sense) a rf-variate Laurent polynomial g, with the same number of terms as / and 
(i = dim Supp(/), such that / vanishes in (R*)" iff g vanishes in (R*)<^ (see Corollary Q of 
Section El below) . In this sense, "almost all" n-variate m-nomials satisfy dimSupp(/) = 
min{n, m — 1}, and those that don't are essentially just rf-variate m-nomials (with d<n) in 
disguise. 

Recall the containments of complexity classes NCi C NC C P C RP C BPP U NP C 
PSPACE (these complexity classes are reviewed briefly in Section|2lbelow). Roughly speak- 
ing, our first main result says that real feasibility for "honest" n-variate (n + /c)-nomials is 
easy for k<2, but NP-hardness kicks in quickly already for k a very slowly growing function 
of n. 

Theorem 1 Let 



A 


:= {feZ[xi,... 




#Supp(/) = 1 + dimSupp(/) andnEN}, 




■■= {/GZ[xi,... 


1 ■^n\ 1 


#Supp(/) = 2 + dimSupp(/)}, 


Ce 


:= |/GZ[xi,.. 




#Supp(/) <n + n,^, dim Supp(/) =n, nGN,l 
and f is a sum of squares of polynomials. J ' 


S 


'■= • • • ) /fc) 


/j GZ[xi, . . . , x„] and fi is a linear trinomial or al 
binomial of degree <2 for alii, and k>n>l. j 



Then, measuring the size of any polynomial f{x) = Yl^=i Cjx"'^ G . . . , Xn] — denoted 
size(/) — as the total number of binary digits in the Ci and Oij, we have: 

1. FEASr(^)gNCi. 

2. For any fixed nGN, FEASM(i3„) gP. 

3. FEASr(5) is NP-hard. 

4. For any fixed 6>0, FEASM(Ce) zs NP-hard. 
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Example 1 A very special case of Assertion (2) of Theorem^ implies that one can decide 
— for any nonzero Ci, . . . , C5 G Z and De'H — whether 

ci + C2X^^'^ + C'iX^'^z + + c^x^y^^z^^ 
has a root in K^, using a number of bit operations polynomial in log{D) + log [(|ci| + !)••• (|c5| + 1)]. 
The best previous results (e.g., via the critical points method, infinitesimals, and rational uni- 
variate reduction, as detailed in IBPROSf ) would yield a bound polynomial in 
D + \og [(|ci| + 1) ■ ■ ■ (Icsl + 1)] instead. Assertion (2) also vastly generalizes an earlier anal- 
ogous result for univariate trinomials \RY05^ . o 

The algorithm underlying Assertion (2) turns out to depend critically on the combinatorics 
of Supp(/), particularly its triangulations. Furthermore, extending the polynomiality of 
FEAS]k(JF) from (n + l)-nomials to (n + 2)-nomials turns out to be surprisingly intricate, 
involving A-discriminants (cf. Section [TTT|) . Baker's Theorem on Linear Forms in Logarithms 
(cf. Section O)), and Viro's Theorem from toric geometry (see, e.g., |GKZ94l Thm. 5.6]). 
Theorem 1, along with some more technical strengthenings, is proved in Section ESI below. 

Remark 1 There appears to have been no earlier explicit statement that FEAS]r(JF) g P 
(or even FEASk(JF) g NP j for T some non-trivial family of n-variate m-nomials with 
m = n-\-0{l). As for lower bounds, the best previous result for sparse polynomials appears to 
have been NP-hardness o/FEAS]r(JF) when T is the family of those n-variate m-nomials 
withm = fl{n^) (see, e.g., IRY051 discussion preceding Thm. 2]). Assertion (4) is therefore 
a considerable sharpening.^ o 

Remark 2 Let Um '■= {/ G Z[xi] | / has exactly m monomial terms}. While it has been 
known since the late 1980's that FEAS^gPSPACE ICanSS^ . it is already unknown whether 
FEAS]r(W4) G BPP U NP, or even whether FEASik(Z[ ]) is NP-hard for some 

particular value ofn ILM01\ [RYO^ . (The latter reference nevertheless states certain analytic 
hypotheses under which it would follow that FEAS^iUm) G P for fixed m.) The role of 
sparsity in complexity bounds for univariate real feasibility is thus already far from trivial, o 

Quantitative results over MJ^ of course shape the kind of algorithms we can find over M". 
In particular, Khovanski's famous Theorem on Real Fewnomials implies an upper bound 
depending only on m and n — independent of the degree — for the number of connected 
components of the real zero set of any n-variate m-nomial |Kho91j . More recently, his 
bound has been improved from 2*^*^™ ^rP^^^rf^"^^ in the smooth case |Kho9H Sec. 3.14, Cor. 



5] to 2°("' )20{n)^o(m) -^^ complete generality |LHWn3[ Cor. 2] (see also |Pern5j for further 
improvements) . 

For n-variate (n + 2)-nomials we can now make a dramatic improvement. Recall that 
a set 5 C is convex iff for any x,y E S, the line segment connecting x and y is also 
contained in S. Recall also that for any AcM", the convex hull oi A — denoted ConvA 
— is the smallest convex set containing A. 



^An earlier version of this paper proved NP-hardness of FEASr(.?") when !F is the family of n-variate 
(6n+6)-nomials, but an anonymous referee suggested an even easier proof for the stronger version of Assertion 
(4) we are now stating. 
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Theorem 2 Let f be any n-variate m-nomial with m<2 + dim Supp(/), Zj^{f) its zero set 
in M" , and define A/'comp(/) (resp. A/'non(/) ) to be the number of compact (resp. non-compact) 
connected components of Z^{f). Then A/'comp(/) < 1 (with examples attaining equality for 
each n&N), and A/'comp(/) = 1 => ^+(/) is either a point, or isotopic to an {n — l)-sphere. 
Also, ATnonl/) is no more than 0, 2, 6, 9, or2n + 2, according as n is 1, 2, 3, 4, or>5 (with 
examples attaining equality for each n<2). Finally, i/Supp(/) is disjoint from the interior 
of its convex hull, then A/'comp(/) =0. 

While the bound Afnonif) ^2 (for n = 2) was found earher by Daniel Perrucci, all the other 
bounds of Theorem |21 are new, and the special case n = 3 improves Perrucci's earlier bound 
for 3-variate 5-nomials by a factor of at least 2 (see |Per05| Thms. 4 & 5]). Except for an 
upper bound of 2'^'^" ^ for the smooth case |Kho9H Sec. 3.14, Cor. 5], there appear to be 
no other earlier explicit bounds in the spirit of Theorem |21 Theorem El has recently been 
generalized to systems of k polynomials in i3„, with identical supports and k E {I, . . . ,n} 
[BRS06J. (See also |BBS05|, iBihOSj for the opposite extreme to Theorem |21 bounding the 
number of isolated real roots of n polynomials in i3„ with the same support.) We prove 
Theorem 121 in Section IT2l 

It thus appears that, unlike algebraic geometry over C, large degree is potentially less of a 
complexity bottleneck over M. Considering the ubiquity of sparse real polynomial systems in 
engineering, algorithmic speed-ups in broader generality via sparsity are thus of the utmost 
interest. Furthermore, in view of the complexity threshold of Theorem Q randomization, 
approximation, and/or average-case speed-ups appear to be the next key steps if we are to 
have a sufficiently general and useful algorithmic fewnomial theory over M. A promising step 
in this direction can be found in work of Barvinok [Barn2] . but more work still needs to be 
done before we can assert significant new randomized algorithms — even for FEAS]k(W4). 

Since the largest m for which FEASk is doable in polynomial-time — for input a single 
n-variate m-nomial — appears to be m = n + 2 (as of early 2006), we propose the following 
conjecture to address the cases m>n + 3. 

Univariate Threshold Conjecture. For any mEN, let 

Um'-= {f ^'^[xi] I / has exactly m monomial terms}. 
Then FEAS]r(Z[xi]) is NP-hard but, for any fixed m, there is a natural probability measure 
onUm so that FEASR(ZYm) has polynomial-time complexity on average. 

Note that FEASM(W3)eP, thanks to Assertion (2) of Theorem^ since U^ = Bi. The latter 
"positive" part of the conjecture is meant to be reminiscent of Smale's 17— Problem, which 
concerns the complexity of approximating complex roots of polynomial systems |SmaOO| 
iRYOHj . 

That feasibility over M may be NP-hard already for univariate sparse polynomials is 
suggested by a recent parallel over a different complete field: Qp. In particular, (a) there is 
now a Theorem on p-adic Fewnomials (due to the middle author j, Roj04| ), with signif- 
icantly sharper bounds than Khovanski's Theorem on Real Fewnomials, and (b) it is now 
known that FEAS(Qp(Z[2;i]) — the natural p-adic analogue of FEAS]r(Z[xi]) — is doable 
in randomized polynomial-time only if an unlikely containment of complexity classes occurs: 
NPCBPP jPEQSj. 
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Possible alternative evidence for the Univariate Threshold Conjecture can be given via 
the A-discriminant, which is defined in Section im below. To set the stage, first recall the 
following elementary example. 

Example 2 Note that f{xi) := a + bxi + cxf has either or 2 real roots according as the 
discriminant A = 1"^ — Aac is negative or positive. Observe then that the real zero set, WgM.^, 
of A can be identified with the collection of all quadratic polynomials possessing a degenerate 
root, and that A also defines a curve W in the real projective plane P^. Furthermore, W is 
equivalent (under a linear change of variables overQ) to a circle, and thus there are exactly 2 
discriminant chambers. These 2 chambers correspond exactly to those quadratic polynomials 
possessing either or 2 real roots. Finally, note that the support of f is {0, 1, 2} and this set 
admits exactly 2 triangulations with vertices in {0, 1,2}.- They are ' ' o,nd I \ ' ■ o 



1.1 Complexity and Topology of Certain ^4- Discriminants 

The connection between computational complexity (e.g., of deciding membership in semi- 
algebraic sets |DL79j or approximating the roots of univariate polynomials |Sma87j ) and 
the topology of discriminant complements dates back to the late 1970's. Here, we point to 
the possibility of a more refined connection between FEASk and discriminant complements. 
(See also |DRRS05] for further results in this direction.) In particular, our last example was 
a special case of a much more general invariant attached to spaces of sparse multivariate 
polynomials. 

Definition 3 Given any A = {ai, . . . , a^} C of cardinality m, define the set 

V° := |(ci, . . . , Cm) gC™ I /(a;) : = X]j=i Cj-a;"^ has a degenerate root^ in (C*)"|. 
The j4-discriminant is then the unique (up to sign) irreducible polynomial A^gZ[ci, . . . , Cm]\{0} 
whose complex zero set contains V° . (If codimVyi > 1 then we set Aa '■= ^■) For 
convenience, we will usually write Aa(/) in place o/ Aa(ci, . . . , Cm)- Finally, we let Va 
— the A-discriminant variety — denote the zero set of A a in C*" . o 



Example 3 // we take A:- 



then AA(ao, . . . , a^, 6o, • • • , K) is 



exactly the classical Sylvester resultant of the univariate polynomials + ■ ■ ■ + a^x and 
&o + ■ ■ ■ + heX^ . This is a special case of a more general construction which shows that any 
multivariate toric resultant can be obtained as a suitable A- discriminant \GKZ94\ The 
Cayley Trick, Prop. 1.7, pp. 274]- o 

Definition 4 Let Ta ■= {/ e Z[xi, . . . ,a:„] | Supp(/) C A}, and let ADISCVAN (resp. 
ADISCSIGNj denote the problem of deciding whether A A{f) vanishes (resp. determining 
the sign of AaU)) for an input feJ^A- Finally, let ADISCVAN(J^) (resp. ADISCSIGN(J^); 
be the natural restriction o/ ADISCVAN (resp. ADISCSIGNj to input polynomials in 
some family T . o 

^That is, a root ( of / with = • • ■ = ^j^^,; = 0. 
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An intriguing link between FEASk and discriminants is the fact that those A with FEASR(^yi) gP 
currently appear to coincide with those A with ADISCVAN(jF4) G P. This should not 
be too surprising in view of the following fact: If /GJF4 has smooth complex zero set (and a 
similarly mild condition holds for its zero set at infinity), then / lies in some connected com- 
ponent C of (M*)*"^ \ Va (under a natural identification of coefficients of / and coordinates 
of (M*)*^), and any other gEC has real zero set isotopic to that of / (see, e.g., |GKZ94| Ch. 
11, Sec. 5A, Prop. 5.2, pg. 382]). Let us call any such C a discriminant chamber (for A). 

So deciding FEASir(JF4) for a given / is nearly the same as deciding whether / lies in a 
particular union of discriminant chambers, and it is thus natural to suspect that the following 
three situations may be equivalent in some rigourous and useful sense: (a) (M*)*"^ \ Va has 
"few" connected components, (b) FEASk(J^^) is "easy", (c) ADI SCVAN( J^a) is "easy". 

While connections between (a) and (c) are known (see, e.g., |BCSS98( Ch. 16]), the 
best current theorems appear to be too weak to yield any complexity lower bounds of use 
for the Univariate Threshold Conjecture. As for stronger connections between (b) and (c), 
concrete examples arise, for instance, from the following two facts: (1) FEASir(Q) G P, 
where Q:={/gZ[xi, . . . ,x„] | / is homogeneous and quadratic} |Bar93j '^ and (2) for A the 
support of a quadratic polynomial /, Aa(/) is computable in polynomial-time. (The latter 
fact follows easily from an exercise in Cramer's Rule for linear equations and the Newton 
identities |BCSS98| Ch. 15, Pgs. 292-296].) 

A new connection we can assert between the "easiness" of FEASr(JF4) and ADISCVAN(jF^) 
is the case of n-variate {n + 2)-nomials, thanks to Assertion (2) of Theorem ^ above, and 
the first part of our final main result below. 

Theorem 3 Following the notation of TheoremU\ and our last two definitions: 



2. For any fixed n, ADISCSIGN(i3„) G P. 

Note that Theorem |31 improves considerably on what can be done through quantifier elimina- 
tion (e.g., |Can88t IBPROSj ) . because for fixed n these older methods already have complexity 
exponential in our notion of input size. Theorem |21 — proved in Section IXD — turns out to 
be a central tool in the algorithms behind the complexity upper bounds of Theorem Q and 
is the main reason that diophantine approximation enters our scenery. 

In light of the connections between the complexity of FEASk(JF^) and ADISCVAN(JF4), 
we conclude our introduction with another possible piece of evidence in favor of the "nega- 
tive" portion of the Univariate Threshold Conjecture: 

Karpinski-Shparlinski Theorem \KS9f)f ADISCVAN(Z[xi]) is computationally hard 
in the following sense: If ADISCVAN(Z[xi]) C C for some complexity class C, then 
NPCCURP.5 In particular, ADISCVAN(Z[xi]) CBPP =^ NPCBPP. ■ 

''Barvinok actually proved the stronger fact that FEASR(Q5fc) e P for any fixed fc, where QSk is the 
family of polynomial systems of the form (/i, . . . , fk) with fi e Q for all i. 

^The paper |KS99j actually asserts the stronger fact that ADISCVAN(Z[xi]) is NP-hard, but without 
a proof. One of the authors of |KS99j (Igor Shparlinski) has confirmed this oversight, along with the fact 
that it was also observed independently by Erich Kaltofen |Shp06| . 



1. ADISCVAN [j Bn]eP. 
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The containment NPCBPP is widely disbelieved, so it would appear possible that ADISCVAN(Z[xi]) 
is not doable in randomized polynomial-time. 

Our main results are proved in Section IHl after the development of some necessary theory 
in SectionElbelow. A useful elementary result on the real zero sets of n-variate (n+l)-nomials 
is then proved in Section El 

2 Background and Ancillary Results 

Let us first informally review some well-known complexity classes (see, e.g., |Pap95| for a 
complete and rigourous description). 

Remark 3 Throughout this paper, our algorithms will always have a notion of input size 
that is clear from the context, and our underlying computational model will always be the 
classical Turing model JPap9^. Thus, appellations such as "polynomial-time" are to be 
understood as "having bit- complexity polynomial in the underlying input size", and the un- 
derlying polynomial and/or 0-constants depend only on the algorithm, not on the specific 
instance being solved. The same of course applies to "linear-time", "exponential-time", etc. o 

NCi The family of decision problems which can be done within time 0(log InputSize), 
using a number of processors linear in the input size. 

NC The family of decision problems which can be done within time poly-logarithmic in 
the input size, using a number of processors polynomial in the input size. 

P The family of decision problems which can be done within polynomial-time. 

RP The family of decision problems admitting randomized polynomial-time algorithms 
for which a ''Yes'' answer is always correct but a ''No'' answer is wrong with 
probability |. 

BPP The family of decision problems admitting randomized polynomial-time algorithms 
that terminate with an answer that is correct with probability at least^ |. 

NP The family of decision problems where a ''Yes'' answer can be certified within 
polynomial-time. 

PSPACE The family of decision problems solvable within polynomial-time, provided a number 
of processors exponential in the input size is allowed. 

Recall also that even the containment P C PSPACE is still an open problem (as of early 
2006). 

A very useful and simple change of variables is to replace variables by monomials in new 
variables. Please note that in what follows, we will sometimes use real exponents. 



^It is easily shown that we can replace | by any constant strictly greater than ^ and still obtain the same 
family of problems |Pap95| . 
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Definition 5 For any ring R, let R"^^^ denote the set of m x n matrices with entries in 
R. For any M = [niij] e R*^^" and y = {yi, . . . , ?/„), we define the formal expression j/*^ : = 
{yT^'^ ■ ■ 'Vri"'^, ■ ■ ■ ,yT^'" ■ ■ -yrT"'")- We call the substitution x:=?/*^ a monomial change of 
variables. ^4/50, for any z:={zi, . . . ,Zn), we let xz:={xiZi, . . . ,XnZn)- Finally, let GL„(Z) 
denote the set of all matrices in Z"^" with determinant ±1 (the set o/ unimodular matrices), o 

Proposition 1 (See, e.g., ILRW03\. Prop. 2].) For any f/, V"gM"^", we have the formal 
identity {xy)^^ = {x^Y {y^Y . Also, if detU j^O, then the function eu{x):=x^ is an analytic 
automorphism o/M", and preserves smooth points and singular points of zero sets of analytic 
functions. Finally, ?7gGL„(Z) implies that e^"'^(R")=R" and that eu maps distinct open 
orthants o/R" to distinct open orthants o/R". ■ 

Via a simple application of Hermite factorization (see Definition IHl and Lemma Q 
below), we can derive the following corollary which reveals why the dimension related hy- 
potheses of Theorem ^ are mild and necessary. 

Corollary 1 Given any n-variate m-nomial f with d = dim Supp{f) <n, we can find (within 
P) a[/GGL„(Z) such that g{y) := f{y^) is a d-variate m-nomial with dim Supp(/) =(i, and 
g vanishes in {M.*Y iff f vanishes in (R*)". Moreover, there is an absolute constant c such 
that size{g) =0{size{fY) . ■ 

Definition 6 UliS!^ URo9^ Given any MgZ™^", the Hermite factorization of M is an 

identity of the form UM = H where U G GLm(Z) and H = [hij] G Z"^" is nonnegative and 
upper triangular, with all off-diagonal entries smaller than the positive diagonal entry in the 
same column. Finally, the Smith factorization of M is an identity of the form UMV = S 
where t/GGLm(Z), V^gGL„(Z), and 5= [s^-] G Z"^^" is diagonal, with Sj,i|sj+i,j+i for alii, o 

Lemma 1 UJIM \EM^ For any M= [mj.j] gZ"^", the Hermite and Smith factorizations of M 
exist uniquely, and can be computed within 0(?t,^ log'^(nmaxjj |mjj |)) bit operations. 
Furthermore, in the notation of Definition 0, the entries of U , V , S, and H all have bit 
size 0{n^ log^(2n + maxjj |mjj|)). ■ 

To prove Theorems ^ and |21 we will first need some tricks for efficiently deciding when 
polynomials in A have roots in (R*)" or R" . 

Lemma 2 Suppose f is an n-variate {n + l)-nomial with affinely independent support A = 
{oi, . . . , On+i} C R", with Oj = (ctij, . . . , CLnj) for all j . Then 

1. f has a root in R" iff not all the coefficients of f have the same sign. 

2. //AcZ" and all the coefficients of f have the same sign, then f has a root in (R*)" iff 
there are indices i,j,j' with — Oij' odd iff there are indices i,j with Oij — Oi^i odd. 

While this last lemma is ultimately elementary, we were unable to find any similar explicit 
statement in the literature. So we supply a proof in Section |3] Another tool we will need is 
a description of certain real zero sets "at infinity" . 
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Definition 7 Given any compact SdW^ andwGW^, the face of S with inner normal w 

— denoted S"^ — is the set of all x^S that minimize the inner product x ■ w. In particular, 
a facet of S is a face S'^ for some w with dim S"" = dim S — 1. o 

Definition 8 The Newton polytope of f , Newt(/), is the convex hull of the support of f . 
Also, for any w = {wi, . . . , Wn) G M", the initial term function of f{x) = YlaeA '^a^" with 
respect to the weight w is In^(/):= ^ CaX"'. o 

aeA 
w-a minimized 

Theorem 4 \LB.WOS\ Theorem 3 and Lemmata I4 and 15] Let f he any n-variate m-nomial f 
with n-dimensional Newton polytope and define A/'(/) to he the numher of connected 
components of f . Then 

1. If Z+{ln^{f)) IS smooth for all weW\ {O}, then Mnonif) < E A/'(In^(/)). 

w a unit inner facet 
normal of Newt{/) 

2. [Af{ln^{f)) > and Z+(In^(/)) smooth] for some u; G M" \ {O} =^ Z+{f) has a 
non-compact connected component. 

3. [f G ^[xi, . . . , Xn] and w = {0, ... ,0,1)] =^ for e > sufficiently small, Z+{f) has 
no more than 7V"(Inu,(/) — e) non-compact connected components with limit points in 
M^-^ X {0}. 

I A/'non(/)>l ^A/'(hi^(/))>l /orsomeu-GM"\{0}. ■ 

The initial term functions above are also known as initial term polynomials or face 
polynomials when the exponents are integral. 

3 The Proofs of Our Main Results: Theorems 3, 2, and 
1 

We will use [k] in place of {1, ... , k} throughout. Let us start by highlighting two of the most 
important theoretical tools behind our proofs: the binomial formula for circuit discriminants 
and Baker's famous result on approximating linear forms in logarithms. 

First recall that AcM" is a circuit iff A is affinely dependent, but every proper subset 
of A is affinely independent.'' Also, we say that A is a degenerate circuit iff A contains 
a point a and a proper subset B such that a&B, A\a is affinely independent, and _B is a 
circuit. For instance, L-^ . and I . . are respectively a circuit and a degenerate circuit. 

We can now summarize what we need to know about A-discriminants when A is a circuit: 

Lemma 3 Suppose A = {ai, . . . , a„+2} CZ" is a circuit of cardinality n + 2. Then, defining 
size(y4) to he the sum of the sizes of the coordinates of A, and letting Z^[f) (resp. Z^{f)) 
denote the zero set of f in (M*)" (resp. (C*)"j, we have: 

^This terminology comes from matroid theory and has nothing to do with circuits from complexity theory. 
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1 . There is a unique ( up to sign ) vector m = (mi , . . . , mn+2) G Z"'^^ such that (a) the coor- 
dinates of m are all nonzero and have gcd 1, (h) X]?=i — 



fc) V"+^m =0 and fd) \m-\-= Voi(Conv(A\{a,})) . . 



^. r/ie vector m from Assertion (1) can he found in time polynomial in size(y4) 



3. Permuting (if necessary) the aj so that rrij > for all j G [k] and mj < for all 
j E{k + 1, . . . ,n + 2}, we have 

/ k \ / n+2 \ / n+2 \ / k 

A^(ci,...,c„+2)= n<0 n o - n ^^rO n^r 

\i=l / \i=k+l / \i=k+l / \i=l 

4. Aa(4, . . . ,^+2) = for some c' = (c'l, . . . , ^+2) ^ (C*)"+2 Z^{f) contains a 
degenerate point ( (and any such ( satisfies = for all i). In particular, 

c=c' 

Z^{f) has at most one degenerate point per open orthant of (M*)". 

Proof: Assertions (1) and (3) of Lemma El follow immediately from |GKZ94[ Prop. 1.8, Pg. 
274]. 

Assertion (2) follows easily from LemmaQ upon observing that m is merely the generator 
of the integral kernel of a suitable integral matrix. 

To prove Assertion (4), observe that by Assertion (3) and Proposition [U Aa{ci, . . . , c„+2) = 
for some c := (ci, . . . , c„+2) G (C*)'^"'"^ =^ c is a smooth point on a hypersurface defined 
by a real binomial equation. Via |GKZ94| Thm. 1.5, Ch. 1, pp. 16], we then obtain that 

{j2^=i Cjx""^^ has a degeneracy C, and any such ( must satisfy the binomial system stated 
above. (In particular, the coordinates are rational in ci, . . . , c„+2 if A generates Z" as a 
lattice, thanks to Lemma Moreover, by Proposition Q and Lemma ^ we easily derive 
that a real binomial system can have no more than 1 solution per open orthant of (M*)". 

Conversely, (^j^""^^ has a degenerate point =^ Aa(ci, . . . ,0^+2) = 0, by the 

definition of A^. ■ 

Baker's Theorem (Special Case) [Bak77] Suppose co, ... ,cn^1, Qi, . . . ,a]veN, C:=max{4,co, . . . ,C]v}, 
and A := log^ max{4, ai, . . . , On}- Also let A := cq + Ci log(a;i) + ■ ■ ■ + Cn log(aAr). Then 
At^O^ |A|>(CA)-(i6^)'"'"'^i°s^. ■ 



3.1 The Proof of Theorem [31 

Theorem (HI will follow easily from the two algorithms we state immediately below, once we 
prove their correctness and verify their polynomial-time complexity. However, we will first 
need to recall the concept of a gcd-free basis. In essence, a gcd-free basis is nearly as powerful 
as factorization into primes, but is far easier to compute. 

Definition 9 '^BS96, Sec. 8.4] For any subset {ai, . . . , oat} gN, a gcd-free basis is a 

pair of sets {{liVi=i, {'^ij} {i,j)e[N]xM) such that (1) gcd(7i,7j) = 1 for all i ^ j, and (2) 
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Theorem 5 \BS9a Cor. 4.8.2 and Thm. 4.8.7 of Sec. 4.8] Following the notation of 
Definition]^ there is a gcd-free basis for {ai, . . . ,aN}, with rj, size(7j), and size(ejj) each 
polynomial in Yle=i size(a£), for all i and j . Moreover, one can always find such a gcd-free 

basis using just O ( (j2^=i^^'^^{<^i)) ] bit operations. ■ 



Algorithm BinomialVanish 

Input: Integers ai, Pi, Ui,Vi, . . . , a^, Pn, U]\f,VN. 

Output: A true declaration as to whether a^^ ■ ■ ■ = ■ ■ ■ . 

Description: 

0. //nil (sign ai)"^ ^ _^ ^^^(sign pi^^ 2 ^^^^ ^^^^^^ < <^^^^ equal. ' ' 
and stop. Otherwise, replace the ai and Pi by their absolute values. 

1. Construct, via Theorem]^ a gcd-free basis {{li}^=i^ {^ij}{i,j)(^[2N]x[fi\) for ai, . . . ,aN, 
Pi, ... , Pn- 

2. //^^^ ejjMj = ^^^^^j^ eijiii /or a// j G [r]] i/ien onij^wt ' 'They are equal.'' and stop. 

3. Output ''They are not equal.'' and stop. 
Algorithm BinomialSign 

Input: Positive integers ai, Pi,ui, vi, . . . , a^, Pn, un, vn- 
Output: The sign of a^*^ ■ ■ ■ — Pl^ ■ ■ ■ P^ . 
Description: 

1. Let C := \og^'^ max{A, ai. Pi, aN, Pn} , M := max{A,Ui,Vi, . . . ,Un,Vn}, and 
E : = I (CM) ~(32JV)*oo'^ c log c ^ 

2. For alliE [N], let Ai (resp. Bi) be an approximation of log ai (resp. \og Pi) within 
(using, say. Arithmetic- Geometric Mean Iteration \Ber03f ).^ 

3. Output the sign of ^X^ili ^j^j j ~ '^i^ij ^'^^ stop. 

Note that while we can certainly compute A := a"^ ■ ■ ■ a^ using a number of arithmetic 
operations polynomial in s: = size(ui) + ■ ■ ■ + size(tiAr), the bit size of A is already exponential 
in s; hence the need for our last two algorithms. 

Lemma 4 Algorithms BinomialVanish and BinomialSign are both correct. Moreover, Al- 
gorithm BinomialVanish runs in time polynomial in 

/:=^i-,^(size(aj) + size(/3j) + size(Mj) + size(fj)) 
and, if N is fixed. Algorithm BinomialSign runs in time polynomial in I as well. 

Proof of Lemma m That Algorithm BinomialVanish is correct and runs in time polyno- 
mial in / follows directly from Theorem El 

That Algorithm BinomialSign is correct and runs in time polynomial in / (for fixed 
N) follows easily from Baker's Theorem: First, taking logarithms, observe that the sign of 

® Other approximation techniques can be used as weU: It is sufficient to use any algorithm that can find 
the 6 leading bits of log within a number of bit operations polynomial in 6 + logiV. 
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a"^ ■ ■ ■ — Pi'^ ■ ■ ■ (3""^ is the same as the sign of 5* := ^Xlili l^S ~ {Yld=i "^j log A 



Clearly then, 
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< E, so Baker's Theorem tells us that 



Step (3) of Algorithm BinomialSign indeed computes the sign of S. So we have correctness. 

To see that Algorithm BinomialSign runs in time polynomial in / for fixed N , first note 
that the Ai and Bi each require 0(log(M)+log(iV)-logE) =0 ((32A^)^°°^Clo g(C) lo g(CM)) 
bits of accuracy.^ So, via our chosen method for approximating logarithms |Bern3j . we see 

that the complexity of our algorithm is polynomial in (iV(32iV)400^Clog(C) log(CM) + Y.t^ log(aiA)) 

_20iNiogN) size(aj) + size(/3j) j Ei=i(size('Ui) + size(fi)), and we are done. ■ 

The Proof of Theorem |3 First note that if our input A is not a circuit, then A is a 
degenerate circuit, and A^(/) is then identically 1. This is because Z^{f) is smooth when 
the support of / is a degenerate circuit (see, e.g., the proof of Case 2 of TheoremElin Section 
I3.2|l . So we can assume that A = {ai, . . . , an+2} is a circuit, / has support A, and that Cj is 
the coefficient of x"'^ in / for all j. 

Thanks to Lemma El Assertion (1) (resp. Assertion (2)) follows straightforwardly from 
the complexity bound for Algorithm BinomialVanish (resp. BinomialSign) we just proved 
in LemmaEJ In particular, the latter lemma tells us that the complexity of ADISCVAN, for 
an input {A, Ci, . . . , 0^+2), is polynomial in Y^^=i log(cjmj) (following the notation of Lemma 
EJ; and the same is true for ADISCSIGN provided n is fixed. The classical Hadamard 
matrix inequality |Mig92| tells us that size(mj) =0(n log(?7, maxj ^{0^^})), so the complexity 
of ADISCVAN is indeed polynomial in size(/); and the same holds for ADISCSIGN 
when n is fixed. ■ 



3.2 Deforming to Polyhedra: The Proof of Theorem [2] 

Let A = {ai, . . . , am} be the support of / and and write f{x) = X]j=i Cjx"^ . Since Z^{f) is 
unaffected if / is replaced by a monomial multiple of /, we can clearly assume without loss 
of generality that ai = 0. 

If n = 1 then / must be a univariate trinomial and Theorem |21 follows immediately from 
Descartes' Rule. So let us assume henceforth that n>2. 



Case 1: A afRnely independent: Letting A' denote the matrix whose columns are 02, . . . , 0^, 

observe that f{y) := f{y^' ^) = Ci + C2I/1 + ■ ■ ■ + Cmym, with m < n. Moreover, thanks 
to Proposition ^ ^+(/) (the intersection of a hyperplane with the positive orthant of M") 
and Z+{f) are diffeomorphic, so we are done. In particular, we see that Z^{f) is either 
empty or a connected, open, C°°, real {n — l)-manifold. ■ 

Remark 4 Note that we can in fact allow arbitrary real exponents for f in our proof above, o 



^Note also that the true number of bits of accuracy we would use in practice is 2'^, where /x is the smallest 
integer with, say, 2'' > (32iV)'*°°^Clog(C) log(CM). This is because while it may be non-trivial to compute 
E exactly, finding 2^ is easy via the old trick of recursive squaring. 

^°One can in fact check in polynomial-time whether ^ is a circuit. See, e.g., Proposition [3 of Section 
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Case 2: yl is a degenerate circuit: Suppose, without loss of generality, that B = {0, 02, ... , a^} 
(with i <m<n + 2) is a circuit, and a£ in the relative interior of Convi? if B intersects the 
relative interior of Convi?. 

Let A" be the n x (m — 2) matrix whose columns are Om-i, ■ ■ ■ ,0,2 and (via Lemma 
define U to be any n x n unimodular matrix U such that UA" is lower triangular, with 
nonnegative diagonal. Defining f{y) '■ = f{y^), it is then easily checked that f{y) is of the form 

Ci + C2y< + ■ ■ ■ + C„l/< with G {0}«+l->^^^{jVm-2} X ^max{j,m-2}-l all j G {2, . . . , m - 1}, 

and G {0}"+2-™ x (Z \_{0}) x Z"^-^. Moreover, Proposition □ tells us that Z+{f) and 
Z+{f) X M"^™" (when Z+{f) is considered as a subset of M^) are diffeomorphic. 

Defining « ....,<.):=< for all J and /i(?/) = -(ci+^2y;2+-+c^-i/^-i) ^ observe then 

-1/a" 

that /iGZ[?/i, . . . ,2/m-i] and, for all y = {yi, . . .,ym)eW^, we have f{y) = fi '"'™=2/m- 
So we see that Z^{f) is exactly the positive part of the graph of the C°° function J^'^'"'" ; 
M™^^ — > M. Clearly then, Z+{f) consists of a union of connected, open, C°°, real (m — 1)- 
manifolds, and Zj^{f)K, Zj^{f) x M"""* is thus smooth and has no compact connected com- 
ponents. So let us now bound Afnonif)- 

If /i is always positive on M™"^ then Z^(f) clearly consists of a single non-compact 
connected component. So we can henceforth assume that Z^{f) is non-empty, which in 
turn implies that every connected component C of Z^{f) has a limit point in R™^^ x {0}. 
Assertion (3) of Theorem HI then tells us that A/'non(/) < A/'(/i — e) for some e>0. So, by 
induction (reducing to lower- dimensional instances of Cases 2 or 3), we are done. ■ 

At this point, we must recall a result of Viro on the classification of certain real algebraic 
hypersurfaces. In what follows, we liberally paraphrase from GKZ94, Thm. 5.6]. 

Definition 10 Given any finite point set A C MJ^, let us call any function uj : A — > M a 
lifting, denote by vr : W^'^^ — > M" the natural projection which forgets the last coordinate, 
and let A:={{a,uj{a)) \ a&A}. We then say that the polyhedral subdivision of A defined 
by {tt{Q) I Q a lower^^ face of Conv A of dimension dim A} is induced by the lifting u, 
and we call S^^ a triangulation induced by a lifting iff every cell of is a simplex. 
Finally, given any f{x) = XlaeA CaX"" G Z[xi, . . . , x„], we define fuj,e{x) := XlaeA CaS'^^'^h"- to 
be the toric perturbation of / (corresponding to the lifting uj) . o 

Definition 11 Following the notation above, suppose dim A = n and A is equipped with a 
triangulation S induced by a lifting and a function s : A — > {±} which we will call a 
distribution of signs for A. We then define a locally piece-wise linear manifold — the 
Viro diagram V(S, s) — in the following local manner: For any n-cell C G E, let Lc be 
the convex hull of the set of midpoints of edges of C with vertices of opposite sign, and then 
define V(T,, s) := [j Lc- o 

C an n-cell 

Example 4 The following figure illustrates 6 circuits of cardinality 4, each equipped with a 
triangulation induced by a lifting, and a distribution of signs. The corresponding (possibly 
empty) Viro diagrams are drawn in the lightest color visible (yellow on the color version of 
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A lower face is simply a face which has an inner normal with positive last coordinate. 
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this paper), o 




Theorem 6 Suppose f{x) — ^^^^4 CgX" ] with Supp(/) = A and dim. A = n, uj 

is any lifting of A, and define Sf{a) =sign(ca) for all a&A. Then for any sufficiently small 
e>0, Z^{f^^^) is isotopic to V(Ej^, s/) \ 9ConvA. In particular, V(T,^,Sf) is a disjoint finite 
union of piece- wise linear manifolds, each possibly having a non-empty boundary. ■ 

Lemma 5 Suppose A is a circuit, is a triangulation of A, n = dim A, and s is any 
distribution of signs on A. Call any point a in the relative interior of A with s{a) opposite 
s{a') for all a' eA \ {a} a caged alternation of {A, s). Then 

1. Z^{f) smooth =^ Z^{f) is isotopic to V(S, Sf) \ ConvA for some S. 

2. V(S, s) has no boundary iff V(S, s) is the boundary of an n-simplex iff [A has a caged 
alternation a and S is the triangulation obtained by the lifting that sends a 1-^ and 
a' for all a' E A \ {a}]. 

Proof of Lemma O By Lemma El and Lemma d it easily follows that A has at most 2 
discriminant chambers in M""*"^, and each such chamber contains a unique toric perturbation. 
Since the topology of Z^{f) is constant on any discriminant chamber containing / (e.g., 
(;KZ941 Ch. 11, Sec. 5A, Prop. 5.2, pg. 382]), we obtain Assertion (1). 

Now note that by definition, any triangulation of a circuit A obtained by lifting has at 
most n + 1 top-dimensional cells (since it is the projected lower hull of a (n + l)-simplex). 
So then, if V(E, s) has no boundary, its convex hull must clearly have dimension n, in which 
case we see that V(E, s) is a union of at least n + 1 simplices of dimension n — 1. So V(E, s) 
is the union of exactly n + 1 simplices of dimension n — 1, equal to the boundary of its convex 
hull, i.e., V(S, s) is the boundary of a n-simplex. This proves the first rightward implication 
of Assertion (2), and the converse is obvious. 

Now if V(S, s) is the boundary of a n-simplex, then A clearly intersects its relative 
interior, which in turn implies that S must be the specified triangulation. (Indeed, it is 
a standard fact that any circuit has exactly 2 triangulations. So the only other possible 
S has exactly 1 n-cell and could not possibly give the V(S, s) we desire.) Furthermore, A 
must then clearly contain a caged alternation, for otherwise V(S, s) would no longer be the 
boundary of an n-simplex. This proves the second rightward implication of Assertion (2), 
and the converse is obvious. ■ 

We are now ready to return to our proof of Theorem El and finish the remaining special case. 
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Case 3: A is a circuit: First note that by Corollary ^ we can apply a monomial change 
of variables and assume ■m = n + 2 without loss of generality. (Moreover, if m<n + 2 initially, 
then every connected component of Z^{f) must be non-compact, thanks to Proposition^) 
Let us first bound Mnonif)'- First note that since A is a circuit, every facet of Conv(A) 
has affinely independent vertex set. So by Case 1, the facet functions of / each automatically 
have smooth zero sets which are either empty or consist of a single non-compact connected 
component. The Upper Bound Theorem of polyhedral combinatorics |Ede87j then implies 
that Conv(y4) has no more facets than a moment n-polytope with n + 2 vertices. The latter 
polytope has exactly 4, 6, or 9 facets, according as n is 2, 3, or 4 |Ede87j . So Assert ion 
(1) of Theorem m then directly implies our stated bound on Mnonif), for n e {3,4}. That 
A/'non(/) < 2 for 72 = 2 foUows from earlier work of Daniel Perrucci using a different argument, 
involving a detailed analysis of the central special case 1 + x + y + Ax°'y^ jPerOSt Thm. 4, 
Assertion (4)]. 

Note in particular that when n= 1, Z^{f) has no non-compact components unless / is 
identically zero. Note also that f{x, y) := {x — l){y — 1) =xy — x — y + 1 has support a circuit 
and exactly two connected components (each non-compact) for its zero set in Rl. So our 
stated bounds for Afnonif) are indeed tight for riE {1, 2}. 

Assume now that n > 5. If Z^{f) is smooth, then Assertion (1) of Lemma El tells us 
that Z+{f) is isotopic to some Viro diagram. As observed within the proof of Lemma EJ a 
triangulation S of A contains no more than d+1 d-simplices, and thus a Viro diagram of the 
form V(S, s) conists of no more than d+1 (rf — l)-simplices. Since any such {d — l)-simplex 
can belong to at most one connected component of V(S, s), and since V(S, s) is a union of 
piece-wise linear manifolds with boundary, we see that V(S, s) \ ConvA has at most d + 1 
non-compact connected components. So by Theorem IHl A/'non(/) <n + 1 in the smooth case. 

Now recall the following two standard inequalities: 

(Nc) KomM < A4omp(/ ~S)+ KomM + 

(A^n) Afnonif) <Kon{f - s) + Afr^on{f + s) 

for e > sufficiently small. Moreover, Zj^{f + e) is smooth for all e > sufficiently small. 
(See, e.g., |Bas99| Lemma 2].) 

Having proved J^nonif) <n + l in the smooth case. Inequality {N„) then immediately 
imphes that Afnonif) <2n + 2, even in the presence of singularities for Z^{f). 

To bound A/'comp(/); let n be arbitrary once again. Observe then that Lemma El implies 
that A/'comp(/) < 1, with equality only if Z+{f) is isotopic to an {n — l)-sphere, as long as 
Z^{f) is smooth. So our bound for A/'comp(/) holds if Z+{f) is smooth. 

On the other hand, if Z^{f) has a singularity, consider first the special case where {A, Sf) 
does not have a caged alternation. LemmaEland (Nc) then imply that A/'comp(/) =0. So we 
have in fact proved a strengthening of the final assertion of Theorem |21 

As for the case where {A, Sf) has a caged alternation, assume without loss of generality 
that it is a„+2- Lemma El then tells us that A/'comp(/) < 1 (with Z^{f) isotopic to an (n — 1)- 
sphere if A/'comp(/) = 1), provided Z^{f) is smooth. So we have our assertion for compact 
components of Z^{f) in the smooth case. 

To conclude, assume that Z^{f) has a singularity. Lemma El then tells us that this 
singularity is unique and that (nr=/ ^T') = jji-mn+2 Y[^^^ So / — e and f + e 

thus lie in opposite discriminant chambers for any £ > 0. (Note that we are still assuming 
that {A,Sf) has a caged alternation.) So one of Z+{f — e) or Z+{f + e) is empty. By Lemma 
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El and inequalities (Nc) and (A^„), we thus obtain that A/'comp(/) < 1 and A/'non(/) = 0. In 
particular, Z^{f) must have consisted of a single point, for otherwise, both Z^{f — e) and 
Z(/ + e) would have been non-empty (by the Implicit Function Theorem). 

Having proved our upper bound for A/'comp(/), we need only exhibit examples proving 
tightness for each n> 1. For n = l there is the obvious example of (x — 1)^ =a;^ — 2x + 1. As 
for n>2, Theorem ini tells us that it suffices to use e{l + xf"" + ■ ■ ■ + x^") — Xi ■ ■ -Xn, for any 
£:>0 sufficiently small. So we are done. ■ 



3.3 Phase Transitions: The Proof of Theorem [T] 

The complexity lower bounds of Theorem — Assertions (3) and (4) — are the easiest to 
prove, so we start there: 



I The Proof of Assertion (3): Recall that 3CNFSAT is the problem of deciding whether 
a Boolean formula of the form B{X) = Ci{X) A ■ ■ • A Cfc(X) has a satisfying assignment, 
where Cj is of one of the following forms: 

Xi V Xj V Xf^, ~'Xi V Xj V Xk, ~'Xi V ~'Xj V A'fc, ~'Xi V ~'Xj V "iXfc; 
i,j, kE [3A^], and a satisfying assigment consists of an assigment of values from {0, 1} to the 
variables Xi, . . . ,X^m which makes the equality B{X) = 1 true. 3CNFSAT is one of the 
most basic NP-complete problems |GJ79j . In particular, for our purposes, let us measure 
the size of a 3CNFSAT instance such as the one above as N. 

Let us now observe that FEAS]r(A^) is NP-hard, where A4 is the family of all polynomial 
systems of the form (/i, . . . , f^) where, for all i, fi G Z[xi, . . . , Xsat] involves no more than 3 
variables and has degree < 1 with respect to each of them.^^ To see why, first observe that 
any Boolean formula B{X) with no occurence of "A" can be converted into a polynomial 
/b(x) gZ[xi, . . . , X37v] via the following table of substitutions: 

Xi h-* Xi 

-iXi 1 — Xj 

For instance, the formula Xi A -1X3 A Xg becomes 

(xi + (l-X3)-Xi(l-X3))+X8-(xi + (l-X3)-Xi(l-X3))x8 = l-X3+XiX3 + X3X8-XiX3X8. 

To any 3CNFSAT instance as above, we can then associate the polynomial system 

(/Ci - l,---,/Cfc - - Xi), . . . ,X3Ar(l -X3Ar)), 

and it is easily checked that B has a satisfying assignment iff Fb has a real root. (Moreover, 
any root of Fb clearly lies in {0, 1}^^.) Note also that the size of Fb is clearly 0{N), and 
that Fb has at least as many equations as variables. Clearly then, FEASir(A^) G P =^ 
3CNFSATGP, and thus FEASk(A^) is NP-hard. 

To conclude, we need only recall that any system of polynomials F chosen from Z[xi, . . . , x„] 
can always be converted to its Shor Normal Form Sp jSho91j . In particular, F has a 
real root iff Sp has a real root, the size of Sp is linear in the size of F, Sp consists of linear 
trinomials and/or binomials of degree <2, and the number of new equations introduced is 
the same as the number of new variables introduced. (More concretely, substitutions like 
X2 — xJ I— ^ (1/1 — xl,y2 — yf,y3 — I/2I/1, ^2 — 1/3X1) can be used to reduce all powers to 2 or less, 

-'^-'^While this construction is well-known in some circles, we have included a detailed explanation for the 
convenience of the reader. 
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and substitutions like yk = CiX"'^ + Cjx"-^ can be used to reduce any polynomial to a collection 
of polynomials, each with 3 or fewer monomial terms.) So the number of variables of Sp is 
bounded above by the number of equations of Sp, and we are done. ■ 



The Proof of Assertion (4): By our proof of Assertion (3), and replacing any system 

(/i) ■ ■ ■ ^ fkj^S with the polynomial /f + ■ ■ ■ + /I, it is easy to see that FEASk(£^) is NP- 
hard, where £ is the family of n-variate lln-nomials that are sums of squares. Indeed, by 
our earlier reduction to 3CNFSAT, we can additionally assume that any polynomial in 
£^ n . . . , Xn] contains — Xi)^, . . . , — as summands. The latter assumption 
easily implies that dimSupp(/) =n for any n-variate polynomial f for then Newt(/) 
must contain a line segment parallel to each and every coordinate axis of M". 

To reduce to the family Cg, simply note that from any f &S and e > 0, we can form the new 
polynomial ge{x, y) '■= f{x) +yf + ■ ■ ■ + y%, where N:= [lOn^/'^] . Observe then that g involves 
exactly n + N variables, g has exactly lln + N monomial terms, and dim Supp((7) = n + N. 
In particular, 

lln + N - {n + N) = 10n< \10n'^/^Y <{n + N)' , 
for all n > 1. Moreover, for any fixed e>0, the size of g is clearly polynomial in the size of 
/. So we then clearly obtain that FEASir(C£) gP =^ FEASk(£^) gP, and we are done. ■ 

The complexity upper bounds of Theorem ^ — Assertions (1) and (2) — then follow 
easily from the two final algorithms we state below. In what follows, we let Z^{f) denote 
the real zero set of / and y4:=Supp(/). 

Algorithm DimPluslCase 

Input: A polynomial f eIjIxi, ... ,Xn] with ajfinely independent support A = {ai, ... ,am} of 
cardinality m, and Cj the coefficient of x"^ in f for all j . 

Output: True declarations as to whether Z^{f) , Z^{f), and Z^{f) are, respectively, empty 
or not. 

Description: 

1. If the Ci are not all of the same sign, then output 

"Z4f), Z^if), and are all non-empty." 

and stop. 

2. If any one of a2 — ai, . . . , — cli has an odd coordinate, then output 
' ' Z^{f) and Z]^[f) are non-empty, but is empty.'' 
and stop. 

3. Output "Z+{f) and Z^{f) are empty.". 

4- IfO^A then output ' ' Z^{f) is non-empty. ' ' and stop. 
5. Output ' ' Z^lf) is empty." and stop. ■ 

Remark 5 Note that dimSupp(/) can be computed in NCi as follows: compute the rank 
of the matrix whose columns are a2 — cti, . . . , — cti via the parallel algorithm of Csanky 
lCsa7(^ . So checking whether a given f is a valid input to Algorithm DimPluslCase can be 
done within NCi. o 
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The Proof of Assertion (1): Here we simply apply Algorithm DimPluslCase. Assuming 
Algorithm DimPluslCase is correct, the complexity upper bound then follows trivially, since 
one can check all the necessary signs and/or parities in logarithmic time, using a number of 
processors linear in size(/). (Note also that we can check within NCi whether an input / has 
Supp(/) affinely independent, thanks to Remark El) So we need only show that Algorithm 
DimPluslCase is correct. 

To prove the latter, note that Lemma |21 implies that the existence of a pair of co- 
efficients of / with opposite sign is the same as Z^{f) being non-empty. So Step (1) 
is correct, and we can assume henceforth that all the coefficients of / have the same 
sign. By Lemma El again, the existence of indices i,j with Ojj — a^^i odd implies that 
Z^{f) \ Z^{f) is non-empty. So Step (2) is correct, assuming its hypothesis is true. 

So let us now assume that has the same parity as Oj^i for all Note that if / 
has a root in then this root must lie in some coordinate subspace L of minimal positive 
dimension. So, by our initial hypotheses, on L, the polynomial / will restrict to an n'-variate 
m'-nomial with m' <n' + 1 and affinely independent support a subset of A. In particular, 
since we now have that all the coefficients of / have the same sign, and since has the 
same parity as Oj^i for all i,j, one final application of Lemma |21 implies that Z^{f) is empty. 
So Step (2) is also correct when its hypothesis is false. 

Steps (3)-(4) are clearly correct since / vanishes at O iff / does not have a nonzero 
constant term. ■ 

Remark 6 Let FEAS^ and FEASj^ denote the obvious analogues of FEASk where we 
respectively restrict to roots in M" or (R*)". Also, paralleling our earlier notation, let 
FEAS]^(jF) and FEASi^(jF) be the corresponding natural restrictions of FEAS^ andFEASj^ 
to inputs in some family T . Our preceding proof then clearly implies that FEAS^{A) G NCi 
(even if real exponents are allowed) and FEASi^(^) gNCi. o 

We can now describe our final algorithm: 

Algorithm DimPlus2Case 

Input: A polynomial f G Z[a;i, . . . , x„] with support A = {ai, . . . , a„+2} of cardinality n + 2, 
dimA = n, ai = O, and an £ G {3, . . . , n + 2} such that B := {oi, . . . , a^} is a circuit and 
lies in the relative interior ofCoiwB if B intersects the relative interior ofConvB. 
Output: True declarations as to whether Z^{f) has compact and/or non-compact connected 
components, along with notification as to whether Z^{f) is a point. 
Description: 

0. Let W denote the set of all inner facet normals, with integer coordinates having no 
common factor, of A. Also let Cj be the coefficient of x'^^ in f for all j , and let A" be 
the n X n matrix whose columns are a„+i, . . . , 02. 

1. If n>2, decide (using Algorithm DimPluslCase j whether there is a w & W with A^ 
affinely independent and Z+(In^(/)) non-empty. If so, then output 

' ' Z^(f) is non-empty and all its connected components are non-compact.'' 
and stop. 

2. If A is a degenerate circuit then do the following: 
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(a) If B intersects the relative interior of CoiwB then do the following: 

i. Find (via Lemma[I\ of Section\^ a unimodular matrix U such that UA" is 
lower triangular and has a nonnegative diagonal with exactly one zero entry. 
Then, replacing f by —f if necessary, assume that all the Cj (except possibly 
Ci) are positive. Finally, define h{x) :=/(x) — Q+ia;"*+i _ . . . _ Cn+2X°'"^^ and 
g{y):=h{y^). 

a. Decide, via a lower- dimensional instance of Algorithm 'DimPlus2Ca.se, whether 
Z^{g) contains at least 2 points. If so, then then output 

is non-empty, smooth, and all its connected components are non- compact. ' ' 
and stop. 

3. If A intersects the interior of ComrA then let Sf{a) :=sign(ca) for all a&A, order the 
sequence of signs Sf{a) in increasing order of a (ifn = l), and do the following: 

(a) Ifn = l and Sf has exactly one sign alternation then output 
' ' Z^{f) is a point of multiplicity 1.'' and stop. 

(b) If {A,Sf) has a caged alternation then do the following: 

i. Decide, via Algorithm Binomia.lYa.nish, whether A a{ci, ... ,Cn+2) =0- If so, 
then output 

' ' Z^{f) has exactly one connected component, and it is a singular point.'' 
and stop. 

a. Decide, via Algorithm BinomialSign., whether sign{A{g)) = (—1)"^"+^ . If so, 
then output 

' ' Z+{f) has exactly one connected component, and it is smooth 
and isotopic to an {n — l)-sphere. ' ' and stop. 

4. Output ' ' Z+{f) is empty.'' and stop. ■ 

It will be useful to observe that the input hypotheses to the preceding algorithm can be 
checked within P: 

Proposition 2 Given any finite set A G 1^"^ with ^A < dim A + 2, we can decide if A 
contains a circuit B, and find the unique such B should it exist, in time polynomial in 
size{A). Moreover, if B contains a circuit and B intersects the relative interior ofConvB, 
then we can find the unique point of this intersection also in time polynomial in size (A). 

Proof of Proposition [2J First, one simply checks via the method of Remark |31 if j^A = 
dimy4 + 2 (for if not, A can not contain a circuit). Then, via Cramer's Rule and the Newton 
identities |BCSS98l Ch. 15, Pgs. 292-296], one simply checks which subsets of A of cardinality 
dim A form facets of A. (Overlaps can be distinguished via a computation of the underlying 
inner facet normals.) If all the facets of A have cardinality dim A, then A itself is a circuit 
and we are done. Otherwise, some facet S of A has cardinahty > dim A, and we use the 
same method recursively to find the unique circuit of S. Via |BCSS98l Prop. 21, Ch. 15, 
pp. 295], it is then easily checked that the bit complexity of this method is no worse than 
0{size{A)^), where size(A) is the sum of the sizes of the coordinates of the points of A. 
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To efficiently check whether B intersects the relative interior of Conv5, we can employ 
any linear programming algorithm with polynomial-time bit complexity^^ as follows: express 
— whenever possible — each 6j G-B as a convex linear combination of points in B \ {h}. If 
some point bi can be expressed in this way, then this point is unique, and we can permute 
the entries of A so that ai = 0, ae = bi, and aj = bj for all j G [i] \ {i}. Otherwise, B does not 
intersect its relative interior. ■ 



The Proof of Assertion (2): First note that / has a nonzero constant term iff / does 
not have O as a root, and this can be checked with just 1 bit operation. So we can assume 
henceforth / has a nonzero constant term. Also note that the polynomial obtained 
from / by setting any subset of its variables to lies in Bn' U A for some n' < n. Moreover, 
since we can apply changes of variables like Xj H-i> —Xi in P, and since there are exactly 3" 
sequences of the form {ei, . . . , En) with Ei G {0, ±1} for all i, it suffices at this point to show 
that FEAS^(i3„) gP for fixed n. 

This will be accomplished by Algorithm DimPlus2Case, and via Proposition |21 and a 
monomial change of variables (employing LemmaQ), we can indeed prepare / to be a suitable 
input to this algorithm. ^'^ So we can now assume that / satisfies these input hypotheses. 
Assertion (2) then follows immediately — assuming that Algorithm DimPlus2Case is correct 
and runs in P for fixed n. So let us now prove correctness and analyze the complexity along 
the way. 

Correctness and Complexity of Steps (O)-(l): The correctness of Steps (O)-(l) follows 

71, _l_ 2^ 

immediately from Theorem |3J Note also that W consists of no more than 
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(n + 2)(n + l)/2 normals, and each such normal can be constructed (employing Cramer's 
Rule and the Newton identities |BCSS98l Ch. 15, Pgs. 292-296]) via n — 1 determinants of 
(n — l) x(n — l) matrices, followed by a gcd computation (see, e.g., jBS96j [Ch. 3] for a detailed 
exposition on near optimal gcd algorithms). Since we've already proved that Algorithm 
DimPluslCase runs in P, Steps (O)-(l) clearly run in P, and no further complexity analysis 
is needed if the hypothesis of Step (1) is satisfied. 

If the hypothesis of Step (1) is false, then (applying Lemma |21 again) this means that 
all the Cj with aj a vertex of ConvA have the same sign. So we can assume without loss of 
generality that the cj are all positive and continue to Step (2). ■ 



Correctness and Complexity of Step (2): First let us define hj{x):=f{x) - Yl^=n+3-j '^i-^"'' 
and gj (y) := hj (y^) for all j G [n + 2 — £] . Note that G ]R[xi , . . . , Xn-j] for all j G [n + 2 — . 

By the Case 2 portion of the Proof of Theorem it immediately follows that Z^{f) is 
smooth and diffeomorphic to the positive part of the graph of —gi (as a function on M"^^). 
Clearly then, every connected component of Z^{f) is non-compact. Furthermore, since 
Ci, . . . , Q-i > (and q > as well, if ae is a vertex of ConvA), Lemma |21 and Assertion (2) of 
Theorem m imply that every connected component must have a limit point on M"^^ x {0}. 
Moreover, it follows easily from the Implicit Function Theorem and Assertion (4) of Lemma 
12] that Z^{gi) has at least 2 points =^ ^+(/) is non-empty. So by induction, Z^{f) is 
non-empty iff Z^{gn+i-i) is non-empty (provided i<n+ 1). 

^^See, e.g., BCSS98, Ch. 15] for a nice description of a barrier method employing Newton's method. 
""^■^We should note that if dim Supp(/) < n initially then CoroUary^implies that every connected component 
of Z+{f) will be non-compact. 
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Note also that by construction, Supp(5'j) is a degenerate circuit for all j G [n + 1 — £] 
(provided i<n+l). So, to simplify notation, we can clearly assume without loss of generality 
that i = n + l and that Supp((?) is a circuit. Just as in the last paragraph, we can still assert 
the two implications (a) Z^{g) has at least 2 points =^ Z^{f) is non-empty, and (b) Z+{g) 
empty ^> Z^{f) is empty. So we are left with the special case where Z+{g) is a point. 
Moreover, unless {B, Sg) has a caged alternation, the coefficients of g will all be positive, 
thus making Z^{f) empty (and this will be correctly declared later in Step (4)). So we can 
assume that Z^{g) is a point and (5, Sg) has a caged alternation. 

Define ge{y) {j^iZ^ Ciy°'i^ +C£y"'" . Assertion (3) of LemmaElthen immediately implies 

that for 6,e>0 sufficiently small, sign{AB{g — S)) = sign{AB{ge))- In particular, g — 6 must 
then lie in the same 5-discriminant chamber as gs, and Theorem IHl then implies that the 
graph of —g attains a maximum value of within In other words, Z^{f) is empty, and 

this would be correctly declared later in Step (4). So Step (2) is correct, and its complexity 
is dominated by an instance in Bn' of Algorithm DimPlus2Case, for some n' < n. 

We can thus assume now that A is a circuit and continue to Step (3). ■ 
Steps (3) and (4): The correctness of Step (3-a) follows immediately from Descartes' Rule 
of Signs, combined with the observation that the derivative of f{x) (or x""^ f{l/x)) is nonzero 
on R+. 

If the hypothesis of Step (3-b) is violated, then all the coefficients of / must have the 
same sign and Z_^_{f) must then be empty (and this will be correctly declared later in Step 
(4)). So we may assume that {A, sj) has a caged alternation. Moreover, we may also assume 
without loss of generality that Ci, . . . , c„,+i > and Cn+2 < 0. 

If Z^{f) is smooth then Assertion (1) of Lemma El implies that Z^{f) is isotopic to one 
of two possible Viro diagrams, easily seen to be either empty or the boundary of an n- 
simplex. Moreover, by TheoremlHl Z+{f) is isotopic to the latter diagram iff sign(AA(/)) = 
sign(A^(/e)) for all e > sufficiently small, where feix) := e 

simple calculation from Assertion (3) of LemmaOlthen tells us that sign(A^(/£)) = (— l)''"n+2 
for all e>0. In other words. Step (3-b-ii) is correct. 

So now assume Z^{f) has a singularity (. Then, Assertion (4) of Lemma El implies that 
( is the only singularity of Z^{f). Moreover, since {A,Sf) has a caged alternation, Z^{f) 
must then be exactly {(}, as already proved toward the end of the Case 3 portion of the 
proof of Theorem 121 So Step (3-b-i) is correct. 

To conclude, observe that any case not satisfying the hypotheses of any of our steps 
results in Z^{f) being empty, and this is correctly declared by Step (4). Furthermore, we 
see that the complexity of Steps (3)-(4) is dominated by a single instance of Algorithm 
Binomial Vanish and a single instance of BinomialSign, for input /. ■ 

4 The Proof of Lemma |21 

As before, we can assume without loss of generality that Oi = O and Ci = 1 by dividing by 
a suitable monomial. Furthermore, we can clearly permute the Oj so that Cj > iff i < A;, for 
some k<n + 1. Let A' be the matrix whose columns are 02, ... , a„+i. Then, via the change 
of variables x= z/ {\c2\, ■ ■ ■ , Icn+il)"^ (which clearly preserves the existence of roots of / in 
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any open orthant of (M*)"), we can then clearly assume that f{x) = 1 + x°'^ + ■ ■ ■ + x"'' — 

Our criteria for checking the existence of roots of / in M" or {M.*)"" then clearly reduce 
to checking whether k <n or whether A' has an odd entry. So let us prove that the latter 
conditions correctly characterize the existence of roots of / in MJ^ and (M*)". 

First note that f{x) = for some xG (M*)" iff 

{•k) =a and 1 + ai + ■ ■ ■ + ak — «/c+i — ■ ■ ■ — «„ = 0, 

for some a = (cti, . . . , a„) G (M*)". Assertion (1) then follows almost trivially: Assuming 
xgM" , the equality k = n and Proposition [T] imply that /(x) = 1 + cti H — ■ + a„ > 0, so there 
can be no roots for / in M". Taking the inverse implication, suppose k<n. Then we can set 

/ \ 

k -\- 1 k + 1 

a:= i i, r r to obtain 1 + «! + ■ — h — a^+i — ■ ■ ■ — = 0. So if we 

can solve x"^' =a over M^, we will have found a root in R" for /. Proposition^ tells us that 
we can indeed (since det A' j^O), so we are done. 

We now focus on Assertion (2). Letting y: = x^, note that 

x^' = a^ y' = y''^'^ = {x^'V = a^, 

thanks to Proposition HJ where S = [sij] is an n x n diagonal matrix with Si^i|s2,2| ■ ■ ■ \ sn,n- 
So we'll be able to find a root in (M*)" for / iff 

Cs?) There are a,yE (R*)" with y^ = and 1 + ai + • • • + afc — a^+i — ■ ■ ■ — = 0. 
Let us now separately prove the two directions of the equivalence in Assertion (2): 
(<^=): First note that the hypothesis is invariant under a common translation of Oi, . . . , a„+i. 
So we can assume ai = O (and = n as given), and our hypothesis then translates into 
N C Z" having at least one odd entry. A having at least one odd entry then implies that 
the mod 2 reduction of A has positive (Z/2Z)-rank, and this in turn implies that si^i is 
odd. (Since left and right multiplication by matrices in GL„(Z) preserves (Z/2Z)-rank, 
and Si^i|s2,2| ■ ■ ■ \sn,n^ Since the map ey(x) := x^ is clearly an automorphism of the open 
orthants of (R*)" (provided V G GL,„(Z)), there must then clearly be some open orthant 
(having exactly j positive coordinates) which is mapped bijectively onto R_ x R^:"^ under 

/ \ 
J + 1 J + 1 



ey. So then define a to be any permutation of the vector 



3 n-j 



\ ^ J 

such that sign(a^) = (—1, !,...,!). Clearly then, y^ = a^ has a solution in (R*)" and thus, 
by iP) and our choice of a, / indeed has a root in (R*)". 

(^^): We will prove the contrapositive. Via translation invariance again, in the notation 
above, we see that our hypothesis is equivalent to all the entries of A' being even, and thus 
all the Si^i must be even. We then obtain, via Proposition Q that = has no roots in 
(R*)" unless a G R" . But then a G R" implies that 1 + ai + ■ ■ • + a„ > (since k = n hy 
assumption), so there can be no roots for / in (R*)". ■ 
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